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We report on Monte Carlo studies of the critical behaviour 
of superfiuid 4 He in the presence of quenched disorder with 
long-range fractal correlations. According to the heuristic ar- 
gument by Harris, uncorrelated disorder is irrelevant when the 
specific heat critical exponent a is negative, which is the case 
for the pure 4 He. However, experiments on helium in aero- 
gel have shown that the superfiuid density critical exponent 
£ changes. We hypothesize that this is a cross-over effect due 
to the fractal nature of aerogel. Modelling the aerogel as an 
incipient percolating cluster in 3D and weakening the bonds 
at the fractal sites, we perform XY-model simulations, which 
demonstrate an increase in £ from 0.67 ± 0.005 for the pure 
case to an apparent value of 0.722 ± 0.005 in the presence 
of the fractal disorder, provided that the helium correlation 
length does not exceed the fractal correlation length. 

PACS numbers: 67.40.-w, 75.40.Mg, 75.fO.Nr, 64.60.Ak 



It is generally believed that the superfiuid transition 
(A-point) of pure 4 He belongs to the classical 3D-XY 
model universality class. Near the critical temperature 
T c , the superfiuid density scales as \T — T c \^ , where ( is 
measured [Q to be 0.674 ±0.003. There has been consid- 
erable interest in studying the superfiuid transition for 
helium in a variety of porous media. Q One such system 
is Vycor which is a random glass with a porosity of or- 
der 30%. Remarkably, the critical exponent £ is found to 
be unchanged by the seemingly huge perturbation repre- 
sented by the glass. This however is the result expected 
according to heuristic arguments of Harris Q| and more 
rigorous recent work [Q which shows that weak, uncorre- 
lated randomness is irrelevant at the unperturbed critical 
point provided that the specific heat exponent is nega- 
tive. The exponent [j] a « -0.026±0.004 is indeed nega- 
tive for 4 He, however Narayan and Fisher || have argued 
that since a is only slightly negative, the crossover to 
the pure 3D XY critical regime is almost logarithmically 
slow. Until recently Vycor was the only porous medium 
in which £ was observed to be unchanged. New experi- 
ments however appear to have added a second material, 
porous gold (which also has short-range correlations), to 
the list. |] 

Remarkable results have been obtained for the aerogel 
system. Aerogel is a fractal silica 'dust' with porosities 
of 95-98% or higher. Despite the fact that aerogel is al- 
most entirely empty space, and that it (unlike Vycor) 



has only a tiny effect on the critical temperature, the ex- 
ponent £ is apparently quite distinctly shifted to 
approximately 0.75 (or larger), indicating that the nom- 
inally weak aerogel perturbation is relevant at the 3D 
XY critical point and aerogel (perhaps) produces a new 
universality class. The specific heat data suggests that 
either hyperscaling is violated or the amplitude ratio is 
exceptionally small. ||[l(J In addition to shifting £ and a, 
aerogel dramatically changes the topology of the 3 He- 4 He 
mixture phase diagram. 11 n| 

It is known that long-range correlations can make dis- 
order relevant, even when a is negative. []l3| , [l4[ Weinrib 
and Halpcrin ||f~3|| have demonstrated this for the special 
case of weak, gaussian distributed disorder with long- 
range correlations. Li and Teitel j[f| have looked at a 
model of non-weak (broadly distributed) but uncorrelated 
disorder and find an apparent increase in the exponent 
£ to a value which depends on the exponent in the alge- 
braic divergence of the disorder strength distribution at 
weak disorder. 

Machta [|6| has considered a model of aerogel as a rela- 
tively uniform medium filled with pores on many different 
length scales. This model was motivated by early exper- 
iments which saw two transitions, one at the usual bulk 
temperature and one at a slightly lower temperature. Re- 
cent improvements in aerogel synthesis techniques appear 
however to have eliminated the larger pores and inho- 
mogeneities so that now only a single transition is ob- 
served at a temperature slightly below the bulk value. |J 
Recently, Huang and Meng have examined a mean-field 
theory in a percolating cluster system. [|~7j 

In order to investigate this fascinating subject, we have 
performed extensive Monte Carlo simulations on the 3D 
XY model for three cases: i) no disorder, ii) uncorrelated 
disorder, and hi) fractal disorder. We consider lattice 
sizes up to 24 3 , and use the Wolff algorithm to mini- 
mize the otherwise severe effects of critical slowing down. 
Pq po[ . The model is defined by compact phase variables 
{9} on sites of a simple cubic lattice 
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We measured the disorder-averaged superfiuid density p s 
as a function of temperature and the system size using the 
usual Kubo formula expression. HI Defining one 'sweep' 
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as growing and reorienting a cluster of spins (on the or- 
der the system size when near the critical point) with the 
Wolff algorithm, typical runs involved 5 x 10 4 warm-up 
sweeps, and 2 x 10 5 production sweeps, with measure- 
ments taken every 200 sweeps. Results were averaged 
over typically 100 disorder realizations. 

A finite-size scaling analysis is crucial to the accurate 
determination of the critical temperature and exponents. 
The scaling ansatz assumes that p s has the form []2l|,|l5{] 
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TL- n G(L/£) = TL- n G[(T - T c )L 1/u ] 



(2) 



where Q = (2 — a)/v — 2. Hence, the dimensionless com- 
bination p s LP" /T is scale invariant at the critical point 
T c . If we assume that hyperscaling holds, then we have 
il = d — 2 = 1, otherwise f2 is a priori unknown. In the 
inset of Fig.® we determine T c = 2.156 ± 0.001 by plot- 
ting p s L/T vs. temperature for uncorrelated disorder 



1 + 5K I>5 ; 5K I>5 G [-A, A] and A < 1. (3) 



of relatively large (but bounded) strength A = 0.7. Sim- 
ilar calculations for the pure (disorder-free) case yield 
distinctly larger value T c = 2.203 ±0.001. In the main 
part of Fig.[j], G is plotted as a function of the scaling 
variable (T - T c )L l / u with the value v = 0.667 ± 0.005 
which gives the best data collapse onto a single universal 
scaling curve. Clearly, this unchanged exponent is con- 
sistent with the Harris criterion and the experimental 
observation in Vycor that uncorrelated randomness does 
not change the universality class. |2^] In order to cross 
check the above result, we have also measured the magne- 
tization m and computed the Binder ratio, j2^j24|] which 
is automatically scale invariant at the critical point, 



< m > 



3 < to 2 > 2 



(4) 



This allows one to estimate the critical exponent v with- 
out assuming a value for fl as was necessary in the case of 
the superfluid density. Plotting U4 against the variable 
(T — T^L 1 /", we found good scaling for the same values 
of T c and v obtained from the scaling of p s . 

We turn now to a discussion of fractal disorder. A ma- 
terial with fractal dimension Df has a mass that scales 
with length like M ~ L Di . If A < 3 then (real) objects 
can never be fractal beyond some finite correlation length 
£ because otherwise the density would vanish. Static 
structure factor measurements p5[ indicate that acid- 
catalyzed aerogel has fractal dimension D{ ~ 2.4 — 2.5 
over a wide range of length scales from ~ 6A out to 
roughly £ ~ 600A Base-catalyzed aerogel (typically used 
in the helium experiments) is believed to have a some- 
what lower fractal dimension and a lower range of length 
scales. Various measurements of the 'fracton' vibrational 
properties of aerogels however place a lower bound on the 
correlation length for the connectivity that is at least an 



order of magnitude larger. (26| The connectivity struc- 
ture may be important because the closed vortex loops 
in the helium are presumably attracted to the aerogel 
strands and are thus sensitive to the connectivity. One 
of the central mysteries shown up by the experiments is 
the following. In the critical regime with reduced tem- 
perature t ~ 10~ 5 the correlation length of the helium is 
expected to be larger (~ 1/im) than even the (estimated) 
connectivity length scale. Nevertheless no evidence of a 
cross-over to the uncorrelated disorder regime is evident 
in the full-pore experiments. Q Very recently however, 
Crowell et. al. [ p7| have performed experiments in the 
regime of lower helium densities where even larger cor- 
relation lengths can be obtained. They find evidence of 
a possible cross-over to the uncorrelated disorder regime 
with a lower value of 

We have considered the possibility that the deforma- 
bility of the tenuous aerogel structure is relevant. On 
length scales beyond the fracton correlation length, aero- 
gel acts to sound waves like a relatively homogeneous 
system with a low speed of sound (~ lOOm/s) despite 
its very low mass density, indicating that its compress- 
ibility is nearly 10 6 times that of ordinary glass. The 
aerogel is known to be sufficiently flexible that it modi- 
fies the collective sound mode dispersion. p8|]§|] On scales 
beyond the fracton correlation length, it is reasonable to 
argue that the aerogel density fluctuations 5$ act as, sim- 
ple, uncorrelated local annealed disorder coupling to the 
magnitude of the helium order parameter in a Ginsburg- 
Landau theory with action 



+ [a + 5 <5$]|V| 2 + [A + fc<S$]|V| 4 . 



(5) 



Integrating out 5<fr for small g,h,k produces only irrele- 
vant couplings. For stronger couplings however, the sys- 
tem can, in the right circumstances, be driven to a tricrit- 
ical point, beyond which the transition is first order. This 
is precisely what happens in 3 He- He mixtures where it 
is a good approximation to treat the 3 He impurities as 
annealed disorder. j2^|l]] This confirms the idea that de- 
formability of the aerogel should be irrelevant. It should 
be noted however that it is probably more physically cor- 
rect to include the constraint J d 3 r5Q — which would 
lead to the slow logarithmic case of Fisher renormaliza- 
tion of the critical exponents. [^9| The possibility that 
this may account for the peculiar features of the specific 
heat data should be looked into in more detail. 

We seem to be left only with the possibility of aerogel 
as quenched disorder whose fractal character extends be- 
yond the lower bound set by the fracton cutoff. A variety 
of schemes have been used to model the aerogel structure. 
[I10Jl2|j3|] We have chose n a simple percolation model 
[30 1 for the fractal structure (probably more appropriate 
for acid-catalyzed than base-catalyzed aerogel) . We gen- 
erate a critical percolation backbone on the 3D lattice by 
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randomly occupying lattice sites with probability p = p c 
and keep only the largest connected cluster. In addition, 
we selected only fractal realizations with porosity in a 
narrow window centered on the median value in order 
to reduce the sample-to-sample fluctuations in the disor- 
der strength. We confirmed that these objects had the 
known fractal dimension |3(| Df ~ 2.5. 

In order to be able to do finite size scaling, while avoid- 
ing the scale dependence of the porosity, a single large 
cluster was generated on an Lo = 48 lattice (giving a 
porosity of about 95%) and divided into smaller sub- 
systems of size L — 8, 12, 16, 24. Simulations were per- 
formed for different subsystems with periodic boundary 
conditions and averaged over the subsystems and differ- 
ent fractal realizations. The bond strength K on the 
fractal was arbitrarily reduced from unity to 0.26. The 
inset of Fig. (||) shows p s L/T vs. T, and appears to give 
a clear fixed point with T c estimated to be 2. 1834 ±0.001 
which is closer to the pure critical point than for the 
uncorrelated disorder model, since the porosity of the 
fractal is so high. In the main part of Fig.(||), we plot 
G vs. the scaling variable and find that the (apparent) 
critical exponent £ increases to 0.722 ± 0.005. We have 
confirmed this result with measurements of the Binder 
ratio. Unlike the case of uncorrelated disorder, we ob- 
served a slow drift downward of U4 with system size at 
the previously determined T c . Taking this out by scal- 
ing the data by the factor U±(T C , 0)/J7 4 (T c , 1/L) yields 
essentially perfect data collapse with v — 0.72 ±0.007, as 
shown in Fig.(||). Assuming a violation of hyperscaling 
gives an anomalous dimension to the superfluid density 
p s oc L-( 1+e ^G{(T - T c )L x /% we can place an approxi- 
mate upper bound \9\ < 0.06. 

It is enlightening to compare the present results to a 
model with disorder of lower dimension, namely infinitely 
long columnar defects. J3l| Recent work on this model 
indicates that v± = £j_ ~ 1 is even larger than for the 
present model. This must be the case in order to satisfy 
the rigorous lower bound Q , since this model does truly 
represent a new universality class, and not simply a cross- 
over. The superfluid density measured parallel to the 
columns has an even larger exponent (z = V\\/v±_ ~ 1.07). 

In conclusion, we have argued that the apparent in- 
crease in the superfluid density exponent in aerogel can 
not be due to a true change of universality class but must 
be a cross over effect in the regime where the helium 
correlation length is less than the (apparently large, but 
necessarily finite) correlation length for the disorder. We 
have performed Monte Carlo simulations in this regime 
for a percolation cluster model of fractal disorder and find 
an increase in the effective exponent to C, = 0.722 ±0.005 
which appears to be roughly consistent with experiment. 
However, we see no apparent violation of hyperscaling, 
and attempts to confirm the unusual behavior of the ex- 
perimental specific heat in our model have proved too 
difficult computationally at this time. 



ACKNOWLEDGMENTS 

It is a pleasure to thank G. Ahlers, T. Witten, D.S. 
Fisher, D. Huse, MW. Chan, C. Hanna, L. Kadanoff, 
N. Read, J. Machta, M. Cha, G. Zimanyi, R. Scalettar, 
R. Singh, S. Tcitcl, J. Reppy, and S. Renn for helpful 
discussions. The work at Indiana was supported by the 
NSF through grant No. DMR-9416906. The work at UC- 
Davis was supported by NSF Grant No. DMR 92-06023 
and by the Los Alamos National Laboratory through a 
LACOR grant. 



[10] Melora Larson, Norbert Mulders, and Guenter Ahlers, 
Phys. Rev. Lett. 68, 3896 (1992). 

[11] S. B. Kim, J. Ma, and M. H. W. Chan, Phys. Rev. Lett 
71, 2268 (1993). 

[12] Recent numerical work has considered g-state Potts mod- 
els [K. Uzelac, A. Hasmy and R. Jullien, Phys. Rev. Lett. 
74, 422 (1995)], a density function model [L. Pricaupenko 
and J. Treiner, op. cit., 430] , a density functional model 
and the Blume-Emery-Grimths model [A. Falicov and A. 
N. Berker, op. cit., 426] to study the effect of correlated 
disorder on lst-order phase transitions and the phase sep- 
aration of 3 He- 4 He mixtures. 

[13] Abel Weinrib and B.I. halperin, Phys. Rev. B 27, 413 
(1983). 

[14] G. Grinstein in Fundamental Problems in Stat. Mech. VI , 
edited by E. G. D. Cohen (North Holland, New York, 
1985). 

[15] Ying-Hong Li and S. Teitel, Phys. Rev. B. 41, 11388 
(1990). 

[16] J. Machta, Phys. Rev. Lett. 66, 169 (1991). 
[17] Kerson Huang and Hsin-Fei Meng, Phys. Rev. B 48, 6687 
(1993). 

[18] Ulli Wolff, Phys. Rev. Lett. 62, 361 (1989). 
[19] W. Janke, Phys. Lett. A 148, 306 (1990). 



[1] G. Ahlers in The Pysics of Liquid and Solid Helium, 

edited by K. H. Bennemann and J. B. Ketterson (wily, 

New York, 1976) Pt. I, Chap. 2. 
[2] John D. Reppy, J. Low Temp. Phys. 87 205 (1992); M. 

H. W. Chan, Physica B 169, 135 (1991). 
[3] A.B. Harris, J. Phys. C 7, 1671 (1974). 
[4] J. T. Chayes, L. Chayes, D. S. Fisher, and T. Spencer, 

Phys. Rev. Lett. 57, 2999 (1986). 
[5] Onuttom Narayan and Daniel S. Fisher, Phys. Rev. B 

42, 7869 (1990). 
[6] Moses Chan (private communication) 
[7] M. H. W. Chan, K. I. Blum, S. Q. Murphy, G. K. S. 

Wong, and J. D. Reppy, Phys. Rev. Lett. 61, 1950 (1988). 
[8] Norbert Mulders, Ravi Mehrotra, Lori S. Goldner, and 

Guenter Ahlers, Phys. Rev. Lett. 67, 695 (1991). 
[9] G. K. S. Wong, P. A. Crowell, H. A. Cho, and J. D. 

Reppy, Phys. Rev. Lett. 65, 2410 (1990); Phys. Rev. B 

48, 3858 (1995). 



3 



[20] The reader is cautioned to note the advice of I. Vat- 
tulainen, T. Ala-Nissila, and K. Kankaala [Phys. Rev. 
Lett. 73, 2513 (1994)], who show that cluster algorithms 
like that of Wolff's severely test the quality of pseudo- 
random number generators. 

[21] Min-Chul Cha, M.P.A. Fisher, S.M. Girvin, Mats Wallin, 
and A. Peter Young, Phys. Rev. B 44, 6883 (1991). 

[22] It is possible that there is a logarithmic correction to 
scaling js| visible in our data, but it is very weak. 

[23] K. Binder, Phys. Rev. Lett. 47, 693 (1981); K. Binder, 
Z. Phys. B 43, 119 (1981). 

[24] We thank D. Huse for pointing out to us that, in order 
to reduce the variance, it is preferable to take the disor- 
der ensemble average of the ratio rather than separately 
ensemble averaging the numerator and denominator. 

[25] Rene Vacher, Thierry Woignier, Jacques Pelous, and Eric 
Courtens, Phys. Rev. B 37, 6500 (1988). 

[26] See for example: D. W. Schaefer, C. J. Brinker, D. 
Richter, B. Farago, and B. Frick, Phys. Rev. Lett. 64, 
2316 (1990), and the many references therein. 

[27] P. A. Crowell, J. D. Reppy, S. Mukherjee, J. Ma, M. H. 
W. Chan, and D. W. Schaefer, Phys. Rev. B (in press, 
1995). 

[28] M. J. McKenna, T. Slawecki, and J. D. Maynard, Phys. 

Rev. Lett. 66, 1878 (1991). 
[29] M. E. Fisher, Phys. Rev. B 176, 257 (1968). 
[30] D. Stauffer, in On Growth and Form, edited by H.E. 

Stanley and N. Ostrowsky (Nijhoff, Dordrecht, 1986), pp. 

79-100. 

[31] Min-Chul Cha and S.M. Girvin, Phys. Rev. B 49, 9794 
(1994). 

FIG. 1. Universal scaling function G vs. the scaling vari- 
able (T - T c )L 1/u for the case of uncorrelated disorder. The 
critical exponent v is estimated to be 0.667 ± 0.005. The in- 
set shows the dimensionless superfluid density p s L/T plotted 
against temperature. 

FIG. 2. Universal scaling function G for the case of 95% 
porosity fractal disorder. The critical exponent v is estimated 
to be 0.722 ±0.005. The inset shows the dimensionless super- 
fluid density p 3 L/T plotted against temperature. 

FIG. 3. Universal scaling curve for the (size-corrected) 
Binder ratio for the case of fractal disorder. The choice of 
the critical exponent v — 0.72 ± 0.007 makes the data for dif- 
ferent system sizes collapse onto a single curve and agrees with 
the conclusion drawn from analysis of the superfluid density. 
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